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Abstract 

Bond breaking in a strong electric field is shown to arise from a 
crossing of the ionic and covalent asymptotes. The specific example 
of hydrogen abstraction from a diamond(lll) surface is studied 
using a cluster model. The addition of nearby atoms in both the 
parallel and perpendicular direction to the electric field are found 
to have an effect. It is also shown that the barrier is not only 
Telated to the position of the ionic and covalent asymptotes. 

I. INTRODUCTION 

Molecular nanotechnology has gained considerable interest in recent years; it’s goal 
is to build complex structures using nanoscale mechanical systems to place (bond) 
reactive molecules with atomic precision. The work of Drexler [1] gives some idea of 
what might be possible when atomic level control of matter is achieved. Positional 
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control has become possible with the development of the scanning tunneling micro- 
scope (STM) and atomic force microscope (AFM). Eigler and Schweizer [2] were able 
to position, with atomic precision, Xe atoms on a Ni surface using an STM tip. How- 
ever, to be truly useful in molecular manufacturing, the atoms must be chemically 
bonded in place and the surface probably has to be more stable and less reactive than 
a metal, as used by Eigler and Schweizer. Therefore it is hard to imagine a first step 
in the construction of any nanocomponent that does not involve removing one of the 
surface atoms. This is clearly true for an etching-like process, but it would also be 
true for adding atoms, where stable surface atoms must be removed to expose a reac- 
tive surface in order to add new atoms. One likely surface atom is hydrogen, which 
is known to tie off the dangling bonds at the surface of carbon or silicon. Removing 
surface atoms can also be important in macroscopic processes as well. For example, 
in diamond growth one of the critical steps is removing the H atoms which terminate 
the growth. Thus there are many reasons to study the selective removal of surface 
atoms. 

Severed years ago Musgrave et al. [3] proposed a tool for abstracting a surface 
hydrogen atom. The disadvantage of this tool it that it must be regenerated after 
each use. Recently Avouris, Nordlander, and co-workers [4] investigated using an 
electric field to remove surface adatoms. This approach seems the most practical 
today for removing surface atoms. Avouris and coworkers identified two mechanisms 
to explain the bond-breaking process: 1) for weak fields the surface-adatom stretching 
modes are excited in steps, via multiple-vibrational excitations, until dissociation 
occurs and 2) for strong electric fields a STM-induced a — » excitation causes 

bond dissociation. The multiple- vibrational excitation process is much slower than 
bond breaking via STM-induced electronic excitation. They also discussed the second 
mechanism in terms of mixing of the a and cr* orbitals resulting in a bond weakening. 
They presented potential energy curves that show that, for the correct orientation of 
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the electric field, the potential becomes increasely shallow, until it has no well and is 
completely repulsive. 

In this manuscript we present an alternative way to view the rapid dissociation 
that occurs in a strong field. We avoid a discussion based on orbitals and only consider 
electronic states, which are physical observables. This approach shows that the field 
induced dissociation can be viewed as a curve crossing caused by the electric field 
preferentially lowering the ionic potential relative to the covalent potential. 

II. METHODS 
A. Cluster models 

The cluster approach is used to study the abstraction of a hydrogen atom from a 
hydrogen passivated diamond(lll) surface. Carbon atoms at the sides and bottom of 
the clusters, which are bonding to more distant carbon atoms in the extended system, 
are bonded to hydrogen atoms. That is, we tie off the dangling bonds with hydrogen 
atoms. The cluster geometries are optimized in the absence of an electric field and 
we do not reoptimize them in the presence of the field. 

To model the effect of the electric field on the bond between hydrogen and the 
diamond(lll) surface we compute potential energy curves for the C-H dissociation 
using quantum chemical methods for five different clusters. The simplest cluster 
contains only one carbon atom, CH4 (see Fig. 1). To this cluster we add the next 
nearest neighbor carbons to obtain C4H10 (see Fig. 1 ). The C4H10 cluster is extended 
vertically by adding an extra double-layer to obtain CioHie (see Fig. 1 ). The C4H10 
cluster is also extended horizontally to include the six nearest neighbor on-top sites 
leading to C13H22 (see Fig. 1 ). Finally the C13H22 cluster is extended vertically by 
adding an extra double-layer to obtain C22H28 (see Fig. 1 ). The smallest cluster, CH4, 
is used as a model to study the effect of using a restricted versus an unrestricted wave 
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function both in the absence and in the presence of a uniform electric field. The larger 
clusters are used to determine the effect of cluster size on the energy barrier for C-H 
dissociation. 


B. Computational details 

All the geometries are fully optimized using the hybrid [6] B3LYP [7] approach, 
in the absence of a field, in conjunction with the 6-31G* basis set [8]. For all the 
clusters, the C-H bond of interest is oriented along the +z axis. The effect of the 
electric field is studied by adding the electric field to the Hamiltonian and obtaining 
a new solution. That is, the orbitals are fully optimized in the presence of the electric 
field. A uniform electric field (f) of strength 0.05 atomic units (1 a.u. = 51.4 V/A) 
is applied, to the optimized structures along the +z axis parallel to the C-H bond. 
We compute the potential energy curves as a function of the C-H distance, in the 
presence and absence of a field, using the 6-31G(d,p), 6-311G(d,p), 6-311++G(d,p), 
6-311++G(2d,2p), 6-31H-+G(2df,2p) and 6-311++G(3df,2p) Pople basis sets [8]. 
All of the B3LYP calculations are performed using Gaussian94 [9]. 

To obtain some insight into the effect of the electric field on bond breaking some 
state- averaged complete-active-space self-consistent-field (SACASSCF) calculations 
are performed of LiH. The Li Is orbital is inactive and the Li 2s and H Is orbitals 
are active. Two 2 £ + states, included in the averaging proceedure, correlate with the 
Li+H and Li + +H~ asymptotes. The 6-31G** basis sets are used. These SACASSCF 
calculations are performed using Molpro96 [10]. 
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III. RESULTS AND DISCUSSION 


A. The LiH model 

The results of the LiH calculations are shown in Fig. 2. The solid lines correspond 
to the results obtained without any electric field. The lower curve dissociates to Li+H, 
while the upper curve dissociates to Li + -fH - . The small hump in the ground state 
potential is an artifact of the simple level of theory, but the curves are qualitatively 
correct. 

When the electric field (f) is added, the energy changes by 
AE=-/if-l/2af 2 -..., 

where /i is the dipole moment and a is the polarizability. The dashed lines are the 
results obtained when an electric field of 0.01 a.u. is added. The ground state is shifted 
to lower energy and the well depth increases. That is, the electric field stabilizes the 
molecule more than the atoms, since the atoms do not have any dipole term. The 
ionic curve undergoes a dramatic change from a bound state potential to an unbound 
potential. The ionic potential crosses the neutral potential at about 12 A. The change 
in the ionic curve arises because the Li + -fH _ dipole moment, which for large distances 
is equal to the separation, and even for a small field, the f/t contribution to the energy 
is very large when the separation between the ions is large. 

The effect of increasing the field to 0.02 a.u. is shown in the figure. We note that 
we have drawn both the adiabatic and approximate diabatic curves to more clearly 
illustrate the ionic/neutral curve crossing. When the field is added, the neutral di- 
abatic potential is shifted to lower energy and the binding energy is increased. For 
bond distances greater than about 2 A, the diabatic ionic potential decreases almost 
linearly with increasing distance, as expected from the /if term. At shorter r val- 
ues, the potential rises rapidly with decreasing bond length due to Pauli repulsion. 
The adiabatic potentials show that the ground state dissociates to the ionic asymp- 
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tote. The barrier for this curve crossing dissociation mechanism decreases as the field 
increases, see the 0.03 a.u. results. Thus if the field is increased to some critical 
strength, the bond will break and the molecule dissociates. 

In our model, dissociation to ions occurs, in this case Li + and H”. It is possible, 
however, that the H ion would loose the extra electron in the strong electric field. It 
is clear that the barrier decreases with the electric field strength, it increases with 
the energy to form the ions (ionization potential (IP) of Li minus the electron affinity 
(EA) of H), and increases with the strength of the bond to be broken. It also depends 
on the differences in polarisability between the molecule and atoms, and between the 
ionic and covalent potentials. 


B. The CH4 model 

We use the CH4 model to investigate the best way to model the electric field en- 
hanced hydrogen abstraction at the B3LYP level of theory. The H3C-H dissociation 
curves are reported in Fig. 3 and are computed using both spin restricted (R) and un- 
restricted (U) B3LYP approaches. At short r the UB3LYP and RB3LYP approaches 
yield the same solution. As the bond length increases, in the absence of the electric 
field, the RB3LYP and UB3LYP curves are very different, since the RB3LYP solution 
dissociates to the ionic limit, CH3~ + H + , while the UB3LYP solution dissociates to 
the ground state systems, CH 3 + H. Clearly a UB3LYP description is necessary to 
obtain dissociation to the correct asymptote in the absence of the electric field. The 
computed D c of 118.5 kcal/mol obtained using a UB3LYP wave function is in qual- 
itative agreement with the accurate value of 112.2 kcal/mol computed by Partridge 
and Bauschlicher [5]. 

In the presence of a strong uniform electric field of strength 0.05 a.u. (2.57 V/A) 
oriented along the +z axis, the H3C-H dissociation does not occur via a radical 
pathway, but via an ionic one. The potential curves show a 1/r dependence for a C-H 
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bond distance greater than 5 A consistent with an ionic description. For a positive 
field and a C-H bond oriented along the +z axis, CH4 dissociates to CH3 + + H while 
for a negative field and a C-H bond oriented along the +z axis, CH4 dissociates to 
CH 3 “ + H+. The CH3 +H + asymptote is much higher in energy than the CH^+H“ 
asyptote because of the large IP of H; as a result, the barrier for C-H dissociation to 
the CH3 +H+ asymptote is much larger. Clearly the H~ asymptote is favored for all 
clusters consider in this work, therefore we restrict our study to this asymptote for 
all larger cluster. 

While LiH near r c and CHj and H“ or CHJ and H+ at infinite separation are 
clearly well described by a closed shell occupation, the question arises, what does the 
wave function look like near the barrier to dissocation. As demonstated in the figure, 
in the presence of a strong electric field, the UB3LYP and RB3LYP potentials are 
very similar, with the UB3LYP potentials being slightly lower in energy. That is, the 
wavefuntion is reasonably well described by a closed shell occupation at all r values. 
Since the deviations between the UB3LYP and RB3LYP potentials are small we use 
a RB3LYP wave function for all the remaining calculations, as it is less demanding 
computationally than the UB3LYP approach. 

The effect of basis set improvement on the barrier is summarized in Table I. 
Expanding the basis set from 6-3lG(d,p) to 6-31lG(d,p) reduces the barrier by 
12 kcal/mol. The population on the dissociating H at the barrier changed from 
—0.64 using the 6-3lG(d,p) basis set to —0.75 using the 6-31lG(d,p) basis set. Thus 
the larger basis set allows more polarization of the charge and stabilizes the ionic 
component of the bonding and therefore lowering the barrier. Addition of diffuse 
functions also reduces the barrier as it is able to describe the distortion in the electric 
field and lowers the ionic potential with respect to the covalent. 
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C. Cluster size effects 


The barrier clearly depends on the size of the cluster used to model the surface, as 
the IP of CH4 is clearly not that of the surface. In addition, the C-H bond strength is 
lowered when spectator H atoms axe replaced by carbon atoms. Therefore we evaluate 
the effect of cluster size by systematically extending the CH4 model in both a parallel 
and a perpendicular direction to the field. The following clusters, as shown in Fig. 1 , 
are considered: C4H10, CujHie, C13H22 and C22H28- 

Before discussing the results, we note a technical problem associated with the 
larger clusters. While the calculations using the basis sets without diffuse functions 
are consistent with our expectations, those with diffuse functions yield Mulhken pop- 
ulation analyses with charges on the dissociating H that can be greater than 2 . It 
appears that in the strong electric field, there is an energy lowering by moving some 
electron density from the cluster to the diffuse orbitals on the hydrogen. These solu- 
tions do not occur for CH4 because of the large IP of CH3. We avoid these solutions 
by not including diffuse functions in the basis set. Since diffuse functions are expected 
to lower the barrier, our computed values should be considered as upper bounds. 

The energy barriers for a field strength of + 0.05 au are summarized in Table I . 
Since the goal of this expansion is to more accurately model the surface, therefore we 
compute some additional properties to see if we can correlate their values with the 
computed barrier. The difference between the ionic and covalent asymptotes should 
be the IP of the cluster without the hydrogen minus the H electron affinity. The 
clusters can polarize in the field, and to account for this, we compute the cluster IP 
in the electric field and the polarisability, a . We also report the C-H bond strength, 
without the electric field, since an accurate model of surface should have the correct 
H- surface bond energy. 

The energy barrier drops by more than 50 % going from CH4 to C4 Hio* This is 
consistent with the weaker C-H bond in C4H10, a lower IP and a higher polarisability 
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for C4H9. Adding an extra double-layer to C4H10 in the direction of the field, to yield 
CioHje, reduces the energy barrier by 50%. As the number of layers increases more 
charge can be polarized towards the hydrogen atom, thus the difference between 
the IP and IP in the field increases. Extending the C4H10 cluster perpendicularly 
to the field, to yield the C13H22 cluster, increases the energy barrier by 30%. The 
C13H22 cluster is not very polarizable in the direction of the field as indicated by 
the small difference in the IP with and without the field. On the basis of the IP 
and of the polarisability, it might have been expected that this cluster would have 
barrier lower than C4H10, but the computed value is larger. The C22H28 cluster 
allows all the interactions described above to occur. The increase in energy barrier 
due to the addition of an extra double-layer is counterbalanced by the decrease due 
to the addition of the surface neighbors. The resulting energy barrier for C22H28 is 
similar to the one for C4H 10 . The difference between the 6-3lG(d,p) and 6-311G(d,p) 
results decreases as the cluster size increases going from 12.6 kcal/mol for CH4 to 5.6 
kcal/mol for C22H28 as basis set requirements are smaller for larger clusters. Increasing 
the number of polarization functions by using the 6-311G(3df,2p) basis set reduces 
the energy barriers by less than 2 kcal/mol. 

The results show that the computed barrier is sensitive to the size of the cluster 
used. The results also show that it is not simply due to changing the IP of the cluster 
and hence changing the location of the ionic and covalent asypmptotes. At the barrier 
the system still have some covalent bonding character and this is not accounted for 
by only considering the asymptotes. The polarisability increases going from CH4 
to C22H28 and does not reproduce the trend in the computed barriers. Our results 
strongly suggest that you must have several nearest neighbors, both in the direction of 
the field and perpendicular to it. On this basis our C22H28 cluster is probably our best 
results, however, this cannot be easily demonstracted because of the computational 
cost of expanding the cluster further. Thus, while it is easy to obtain insight into the 
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mechanism, it is very difficult to compute quantitative barriers. 


IV. CONCLUSIONS 

In the absence of an electric field the C-H dissociation occurs via a radical pathway. 
In the presence of a positive electric field the C-H dissociation occurs via an ionic 
pathway to yield R + + H~. Cluster size and basis set effects play an important role 
in the determination of the energy barriers. 
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TABLES 

TABLE I. The C-H bond energy, the ionization potential of cluster minus a hydrogen, 
the same IP but in a uniform electric field, and energy barriers, for the abstraction of 
hydrogen in a uniform electric field. The used field strength is +0.05 au (2.57 V/A). The 
polarisability, a, is in bohr 3 . All the remaining quantites fire in kcal/mol. 


Cluster 

C-H 

IP 

IP(f) 

a 

Barrier 


6-31G* 

6-31G* 

6-31G* 


6-3lG(d,p) 

6-31lG(d,p) 

CH£ 

119.3 

247.1 

254.9 

12.71 

78.4 

65.8 

C 4 H 10 

107.4 

184.7 

188.0 

40.94 

36.6 

26.0 

CioH 16 

106.4 

164.0 

185.1 

94.86 

18.1 

10.2 

C 13 H 22 

107.3 

172.4 

172.9 

107.48 

50.5 

43.1 

C 22 H 28 

107.5 

155.6 

163.5 

181.15 

35.3 

29.7 


a The barriers in the 6-31l++G(d,p) and 6-311++G(2df,2p) basis sets are 53.6 and 
54.6 kcal/mol. 
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FIGURES 


FIG. 1. The B3LYP/6-31G* optimized cluster geometries. 

FIG. 2. The SACASSCF potential for LiH as a function of electric field strength. For 
the 0.02 au field we have drawn approximate diabatic potentials as well as the adiabatic 
potentials. 

FIG. 3. The B 3 LYP potentials for the dissociation of an H from CH 4 as a function 
of electric field strength. Both the spin restricted (R) and unrestricted (U) potentials are 
given. 


14 






Energy, kcal/mol 


150.0 

100.0 

50.0 

o.o 

- 50.0 

- 100.0 

- 150.0 

- 200.0 







Energy, kcal/mol 



0.9 1.9 2.9 3.9 4.9 5.9 6.9 

H 3 C— H distance, Angstroms 



